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1. Introduction 

A theoretical framework, based on tlie idea of effective field theory, for describing 
nuclear system was proposed at the beginning of 1990's by Weinberg [1]. The framework 
provides us with a systematic scheme for constructing nuclear operators such as nuclear 
forces and nuclear electroweak currents from an effective Lagrangian of the underlying 
theory. We will refer to it as nuclear effective field theory (NEFT). For reviews, see, 
e.^f., Refs. [2, 3, 4, 5, 6] and references therein. Meanwhile, phenomenological models 
for the nuclear operators have been conventionally used in nuclear physics. Expressions 
and behaviors of those nuclear operators are quite different among the phenomenological 
models and NEFT, whereas all of the nuclear operators give essentially the same reaction 
rates for low-energy reactions in few-nucleon system, e.g., solar- neutrino reactions on the 
deuteron[7, 8, 9, 10]. This implies that there exists an equivalence relation between the 
nuclear operators of phenomenological models and NEFT, and thus it is interesting to 
seek for the relation. 

Recently, one of the authors (SXN) argued the relation for the nuclear forces from a 
viewpoint of the renormalization group (RG) [11]. In that work, the author formulated 
a "scenario" , which we describe in detail in the foUowings, for the relation between the 
nuclear forces based on NEFT (Veft) and phenomenological nuclear forces (Vph) paying 
attention to a difference in the size of the model space between Veft and Vph- One 
may construct many different Vph all of which reproduce the low-energy A'"A'"-data. Those 
nuclear forces are different from one another in modeling small scale phenomena. As high- 
momentum states of the nucleon are integrated out and thus the model space of Vph is 
reduced, information relevant to the details of the small scale physics arc gradually lost. As 
a result, all Vph eventually converge to a single nuclear force defined in a sufficiently small 
model space; the model dependence found in the original Vph disappears.^ The short-range 
part of the single low-momentum interaction is accurately parameterized with the use of 
simple contact interactions, which is the way NEFT describes the small scale physics. 
After all, the parameterization of the single low-momentum interaction is nothing more 
than Veft- By construction, Veft obtained in this way does not have a dependence 
on modeling the small scale physics. This is the scenario for the relation among Veft 
and Vph- In Ref. [11], the author employed a Wilsonian renormalization group (WRG) 

unique low-momentum effective A'^A^ potential with the cutoff value A = 2.1 fm ^ , known as Viow-k, 
is derived using the Lee-Suzuki method by Bogner et al.[12]. 
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equation derived by Birse, McGovern and Richardson [13]^ for the model space reduction, 
and showed a result to convincingly argue that the scenario for the relation between Vph 
and Veft is realized. 

It would be interesting to extend the WRG analysis to the study of nuclear electroweak 
currents. This is the main subject of this paper. Similar to the nuclear forces, the 
expression and behavior of a nuclear exchange current operator of a phenomenological 
model are quite different from those of another model or NEFT, and thus the reaction 
mechanism of an electroweak process is also model dependent. The model dependence 
stems from modeling small scale physics. Thus, if one reduces the model spaces of the 
different nuclear current operators, it is expected that all the operators evolve to be 
essentially a single nuclear current, which may be accurately simulated by the NEFT- 
based parameterization. In this way, nuclear electroweak currents from phenomenological 
models and from NEFT are related through RG. 

In this work, we derive a WRG equation of nuclear current operator for two-nucleon 
system^. By using the WRG equation, we reduce the model space of nuclear current 
operators of either phenomenological models or NEFT with the pion, and examine the 
evolution of the operators. We are specifically concerned with the axial vector current 
associated with the Gamow-Teller transition in the low-energy neutral-current neutrino 
reaction on the deuteron {ud J^pn), where our analysis can be simplified due to the 
absence of the Goulomb interaction. Finally, the effective operator, obtained as a result of 
the model-space reduction, is simulated by a NEFT-based operator without the pion [15]. 

Here, we mention a choice of the model-space reduction scheme. So far, some methods 
have been used to derive a low-momentum A^A'^ interaction from phenomenological A^A'^ 
interactions. Those methods are the WRG method[ll, 13], the Lee-Suzuki method[12, 16], 
and the unitary transformation mcthod[17]. For detailed comparison of them, see Ref. [18]. 
In Ref. [11], the author argued that the use of the WRG method is consistent with the 
construction of the effective Lagrangian in NEFT. This argument is based on that an 
effective Lagrangian in NEFT is in principle obtained by integrating out high-energy 
degrees of freedom within the path integral formalism, and that solving the WRG equation 
is equivalent to performing a path integral. Also in the reference, the three model-space 

''^The WRG equation derived by Birse et al. is equivalent to the Bloch- Horowitz method[14]. 

^We will restrict ourselves to two-nucleon system in this work, and thus consider neither the genuine 

three (or more)-body currents nor the generation of them due to the model-space reduction. 
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reduction schemes are applied to the nuclear forces and it is shown that only the WRG 
method generates an effective low-momentum interaction which is accurately simulated 
by the NEFT-based parameterization of the nuclear forces, even in a case of a small model 
space relevant to a pionless EFT [EFT(]f)].^ Therefore, we regard the WRG method as 
being the most appropriate for studying the relation between the models and NEFT, and 
we adopt it as the model-space reduction method in this work. 

This paper is organized as follows. In Sec. 2, we derive the WRG equation of the nuclear 
current operators for two-nucleon system. In Sec. 3, we present the current operators 
employed in this work, and give a detailed description of our numerical analysis. In 
Sec. 4, we show results of the numerical calculation and give discussion on our results. 
Finally, we give a concluding remark in Sec. 5. In appendices, we give a detailed derivation 
of the WRG equation, a discussion of a conservation law for an effective current operator, 
and a discussion on a role of the tensor nuclear force in a pionless theory. 

2 Wilsonian renormalization group equation for nuclear current operator 

We briefly discuss a Wilsonian renormalization group (WRG) equation for nuclear 
current operators. A detailed derivation of the WRG equation is given in Appendix A. 

We start with a matrix element of a current operator in a two-nucleon system evaluated 
in the momentum space. 



nucleon state, which are derived from an effective interaction with a cutoff A (see Eq. (14) 
in Appendix A) . The quantity p {p') is an on-shell relative momentum for the initial (final) 
two-nucleon state, p = y/ME (p' = y/ME') with E (£") being the energy for the relative 
motion of the two nucleons and M is the nucleon mass, and a (/3) specifies a partial wave. 
The radial part of the current operator between the /3 and a partial waves are denoted by 
0'^^'°'\k', k). The quantities k and k' are the magnitudes of the relative off-shell momenta 
of the two-nucleon system. The cutoff A, which is the maximum magnitude of the relative 
momenta, specifies the size of the model space spanned by plane wave states of the two- 
nucleon system. 

^An effective interaction obtained witli the WRG method is on-shell energy dependent. The impor- 
tance of considering the on-shell energy dependence of a NEFT-based nuclear force has been examined 




(1) 



where '^'^\k) {ijjp,{k')) is the radial part of the wave function for the initial (final) two- 



in Ref. [11, 19]. 
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Figure 1: A graphical representation of the WRG equation for a current operator. O (O') 
is an original (renormalized) current operator and is a iViV-interaction. In the loops, 
the magnitude of the relative momentum of the two nucleons is denoted by A;, which lies 
in the momentum shell that is integrated out. 

We differentiate the both sides of Eq. (1) and impose a condition that the matrix 
element is invariant with respect to changes in A, i.e., d{0)/dA = 0. This leads to the 
WRG equation for the effective current operator/ 



where V^" (/c', k]p) is an effective iViV-potential for a partial wave a, which evolves fol- 
lowing the WRG equation for a A^A^ potential [11, 13]. In vl"'\k' , k; p) , k (k/) denotes 
the off-shell relative momentum of the two nucleons before (after) the interaction, p de- 
notes the on-shell momentum, and A is the cutoff value specifying the model space. Some 
arguments in the current operator, which have been suppressed in Eq. (1) for simplicity, 
are shown. We note that, because the current operator and NN potential are evolved 
by means of the WRG equation, the effective current operator acquires a dependence on 
both the initial and final on-shell momenta, p and p' {E and E'), whereas the effective 
A^A^-potential acquires a dependence on either the on-shell momentum of the initial state 
or that of the final state, p or p' {E or E'). A graphical representation of Eq. (2) is 
shown in Fig. 1. For an infinitesimal reduction of the cutoff, the operator O evolves into 

In Ref. [11], the WRG equation for the nuclear force is derived using the path integral method. The 
WRG equation for the current operator, Eq. (2), is also derived in essentially the same way. 



dA 



M ( 0[^^''\k',A;p',p)vt\A,k;p) , vl^\k',A;p')0[^^''\A,k;p',p) 



(2) 
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the renormalized one O' by absorbing the one-loop graph either before or after the in- 
sertion of the current operator. The loops include the intermediate two-nucleon states of 
A — < A; < A, where k is the magnitude of the relative momentum. The WRG equation 
given by Eq. (2) is for a single A^A^ partial wave, but an extension to a coupled- channel 
case is straightforward. 

One starts with O = 0\ for a cutoff A and derives Oa for A(< A) by solving the WRG 
equation. A solution (integral form) of the WRG equation, Eq. (2), may be written as 

Of"^ = f of + of "^A ^-^XVt^ + vfx ^^AOf") 



+ ^f^-^ ^sY^C'f^^A h^Xvt^\r}, (3) 



where H^^\ V^"''^ are the full Hamiltonian and the A^A^-interaction for a partial wave a, 
respectively, defined in the model space with A. Furthermore, r) and A are the projection 
operators defined by 



V 
A 



k<A, (4) 



J^\k){k\, A<k<A, (5) 



where |A;) represents the radial part of the free two-nucleon states with the relative mo- 
mentum k. By inserting Eq. (3) into Eq. (2), one can easily check that the integral form 
of the WRG equation, Eq. (3), satisfies its differential form, Eq. (2). 

A nuclear current operator has the properties which a nuclear force does not possess, 
i.e., the current operator satisfies a (partial) conservation law. It is interesting to examine 
the conservation law for the effective current operator obtained after the model-space 
reduction. This subject is discussed in Appendix B. 

3. Nucleetr current operators 

In this section, we introduce the nuclear current operators used in our numerical renor- 

malization group analysis. We employ five nuclear current operators contributing to the 

GT transition: two of them are phenomenological models and the others are based on 

pionful EFT [EFT(7r)] with three cutoff values ^. Each of the five operators has the same 

^ The EFT(7r)-based operators contain a Gaussian type cutoff function. We denote the cutoff by Ag/t 
and distinguish it from the sharp cutoff A introduced in the WRG equation. 



6 



A 



Jt(p) 



(a) 
/7s 



jt(p) 



(b) 



t(p) 



(c) 



Jt(p) 



t c P 



(d) 



(e) 



Figure 2: Diagrams for axial- vector exchange cur- 
rents of phenomenological models. Upper two di- 
agrams (a,b) are for the A-excitation currents in- 
duced by the tt (or p) exchange. Lower diagrams 
(c) and (d) are for the tt (or p) exchange pair 
currents, and diagram (e) is for the 7r-p current. 
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Figure 3: Diagrams for axial- vector 
exchange currents from EFT(7r). 
Dashed lines denote pions and ver- 
tices without (with) white box arise 
from the leading (sub- leading) order 
Lagrangian. 



one-body operator, 

^i(^) = -9AY.t5i5^'\x-ri). (6) 
i ^ 

where qa is the axial- vector coupling constant, and Tj and iTj are the isospin and spin 
operators, respectively, for the i-th nucleon. The isospin operator has only the third com- 
ponent, because we consider the neutral-current reaction, vd — > vpn, in this work. 
It is noted that we neglect the small corrections from the axial form factor, the induced 
pseudoscalar current, and corrections from higher order terms {e.g., 1/M corrections) be- 
cause we are concerned with the reaction in a low-energy region, Ei, < 20 MeV, where Ei, 
is the incoming neutrino energy. 

We now present the two-body operators. The operators of two phenomenological mod- 
els are found, e.g., in Ref. [8]. The operator for one of the models are obtained by 
calculating five diagrams shown in Fig. 2. The model was originally proposed in Ref. [20] 
and the strength of the A-excitation currents in the model was fixed to reproduce the 
experimental tritium /?-decay rate. We will refer to the model as "Model I" . The other 
model consists of only the A-excitation currents (the upper two diagrams in Fig. 2). The 
overall strength of the model was adjusted so as to reproduce the same total cross section 
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for ud upn as that of Model I at the reaction threshold. We will refer to this model as 
"Model II". For explicit expressions of the operators, see Eqs. (22)-(26) of Ref. [8]. 

The nuclear currents based on EFT(7r) are obtained from heavy-baryon chiral La- 
grangian up to next-to-next-to-next-to leading order in Refs. [10, 21]. Diagrams for the 
two-body current operator are shown in Fig. 3 and explicit expressions for the operators 
are given in, e.g., Eq. (20) in Ref. [21]. We use three EFT(7r)-based operators with dif- 
ferent cutoff values, Ag/t = 500, 600, and 800 MeV (see the footnote 8). For each of the 
cutoff values, the coupling constant of the axial-vector-four-nucleon contact interaction 
shown by diagram (d) in Fig. 3 has been fixed so as to reproduce the tritium /?-decay 
rate [22] 

In our numerical analysis, we consider the operator, O, obtained with the current 
operators, A{x), presented above, 

= J dx e"^-^A{x) , (7) 

where A{x) = yl(a;)|^^Q, with A{x) being the axial-vector current. The center of mass 

— * 

coordinate of the two nucleons is denoted by R. The quantity q is the momentum transfer 

to the nuclear current operator. We will examine the evolution of the operator O defined 

in Eq. (7), rather than the current operator A{x) itself, using the WRG equation. Also, 

a matrix element we calculate is always that of the operator O. In what follows, we refer 

to the operator O in Eq. (7) also as "current operator". 

Matrix elements of the operators are calculated before a model-space reduction (A = 

oo), and the numerical results are shown in Table 1. We calculate contributions (ratio) 

to the matrix elements from the one-body and the two-body current operators. We 

decompose the contributions from the two-body operators into those from the deuteron 

S and D-wave states. For the final neutron-proton scattering state, we consider only 

the ^5*0 partial wave, which is the dominant state in the low-energy region. The total 

amplitude from "Model I" is normalized to 100%. In calculating the matrix elements, we 

employ a specific kinematic condition: E — —B, E' — 1 MeV, and q — ?>Q MeV where 

B is the deuteron binding energy and q — \q\}^ The wave functions are obtained with 

^Recently, it is pointed out in Ref. [23] that a contribution from the contact interaction (Fig. 3(d)) 

plays an important role to deduce a value of neutron-neutron scattering length a„„, from experimental 

data of the 7r~rf nwy reaction. The authors fixed the coupling of the contact interaction in a cutoff 

independent manner, so as to reproduce &pp ^ de+i/g reaction matrix element from a model calculation. 
^•^The total cross section of the vd — > vpn reaction at = 20 MeV gets a good amount of contribution 
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the Argonnc vl8 potential [24]. Since all of the five operators have the same one-body 
operator, we have the same contribution from the one-body operator with the deuteron 
S- and D-waves, 98.87%. Meanwhile, the five two-body currents give small contributions 
which agree with each other within 1% accuracy. We have 1.13, 0.72, 1.42, 1.41, 1.40% 
contributions from the two-body currents with the deuteron S- and D-waves from the 
Model I and II, and EFT(7r) with the three cutoff values A^ft = 500, 600, and 800 MeV, 
respectively. 

We can see from the table that the behaviors of the two-body currents are quite dif- 
ferent. Model I and II induce the "^Di ^So transitions (1.47 and 0.72%) more strongly 
than the ^5*1 ^Sq ones (—0.33 and 0.00%), because the important operator in the 
models is the tensor-type originated from the A-excitation currents exchanging a pion 
or a rho meson. In contrast, the EFT(7r)-operators generate the ^Si — > ^So transition 
amplitudes (-1.18, -2.00, and -3.26% for A^ft = 500, 600, and 800 MeV, respectively) 
and ^Di ^So ones (2.60, 3.41, and 4.65%), both of which are larger than those of Model 
I and II. The ^Si — > ^-S'o transition is due to the contact interaction, the diagram (d) in 
Fig. 3, whereas the ^Di — > ^5*0 transition is due to the pion exchange two-body currents, 
mainly the diagram (c) in Fig. 3. We find that a larger cutoff value (Ag/t) leads to a more 
destructive interference between the deuteron S and D-wave contributions. The interfer- 
ence between the two transitions are so destructive for the EFT(7r)-based operators and 
thus the net contribution is almost the same as those of Model I and II. 

4. WRG analysis of the current operators 

Starting with the five sets of the nuclear operators discussed in the previous section, we 
calculate effective operators at the cutoff values A = 200 and 70 MeV using the integral 
form of the WRG equation, Eq. (3). In this work, we define the one- and two-body 
effective currents as 



01^^ = 7]Oi7], (8) 

L Jii — n 

1 1 1 

V, (9) 



+^A^A(0, + 0.)^AI^; 



where Oi and O2 are the "bare" one- and two-body current operators, respectively, which 
we have already introduced. We note that the high energy parts of the bare one-body 



from the kinematical region around this kinematics. 
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(a) (b) (c) (d) 




Figure 4: A graphical representation of the renormahzation of the effective two-body op- 
erator due to the one-body operator. In the diagrams (a)-(c), all external two-nucleon 
systems have the relative momentum belonging to the model space while all the interme- 
diate states are the high momentum states to be integrated out. After the model-space 
reduction, the diagrams (a)-(c) contribute to the renormalization of the two-body operator 
shown by the diagram (d). 



current operator Oi are included in ^ , as graphically shown in Fig. 4. The cutoffs in 
the projection operators, rj and A in Eqs. (4) and (5), are thus K = oo and A = 200 or 70 
MeV. 

In our WRG analysis, we observe that the effective two-body operators suddenly jump 
up at a momentum around the cutoff. (We will see it below.) As shown in Eq. (9) 
(graphically in Fig. 4), the bare one-body current operator can give a contribution to 
the effective two-body operator, and this is the origin of the "jump-up" structure in the 
effective operators. In Fig. 4(a) (Fig. 4(b)), the high momentum states before (after) 
the insertion of the current operator are integrated out whereas, in Fig. 4(c), the high 
momentum states before and after the insertion of the current are integrated out to 
renormahze the two-body operator. The diagrams (a), (b) and (c) in Fig. 4 correspond to 
the Oi part of the second, third and fourth terms in Eq. (9), respectively. Because of the 
momentum conservation, these terms of Figs. 4(a)- (c) give non- vanishing contributions 
only if the condition, ^ — | <^'<^-|-|, is satisfied; k (k') is the magnitude of 
the relative momentum of the two-nucleon state just before (after) the insertion of the 
current operator. The diagrams. Figs. 4(a) and (b), start generating a large amount of 
the contributions at a certain kinematical point. For example, the Oi part of the second 
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term in Eq. (9) gives noii- vanishing contribution for 

<k' <A. (10) 



A-f 



The behavior of the renormahzed effective two-body operators thus changes suddenly at 
the value of the momentum k {k'), as discussed above. 

4.1 Effective current operators at A = 200 MeV 

Now we are in a position to discuss our result of the effective operators for A = 200 
MeV. We reduce the model space of Model I or II or EFT(7r) using the WRG equation, 
and evaluate the GT matrix elements for the kinematics specified before. The result is 
shown in Table 1. The contribution from the effective one-body operator, 97.05%, is 
smaller than that (98.87%) for A = oo by 1.82%. This is because a part of the one-body 
operator sensitive to the high momentum part of the wave functions is renormahzed into 
the effective two-body operator. In addition, we find that most of the two-body D-wscve 
contributions at A = cxo arc also renormalizcd into the two-body S'-wave ones at A = 200 
MeV. We can see that there still remains some model dependence among the models and 
EFT(7r) at A = 200 MeV. On the other hand, the three EFT(7r)-based operators, which 
appear differently at A = oo, give almost the same result at A = 200 MeV. 

In Figs. 5 and 6, we show the radial part of effective two-body current operators for A = 
200 MeV obtained from the EFT(7r)-bascd operators. We denote the effective operator 
for the initial deuteron S- and D-wave states (and the final ^Sq state) by 02^s{^' ^ ^) 
and 02^g^{k' , k), respectively. Wc plot the diagonal part of them, i.e., k' = k, and also 
that of the "bare" current operators, 02,ss{k, k) and 02,sD{k, k). Wc find that the curves 
of the effective operators in the figures suddenly change at k = A — q/2 = 185 MeV 
("jump-up" structure). This is because, as discussed above, at this point the diagrams 
Figs 4(a) and (b) start contributing to the effective two-body operators. We also find 
that the cutoff Aeft dependence among the original EFT(7r)-based operators at A = oo 
are much less observed at this resolution of the system. However, the disappearance of 
the model dependence is not as perfect as the case for the nuclear force, Viow-k, in which 
the model dependence of phenomenological nuclear forces are almost not observed at 
A(= 2.1fm~-^) ~ 400 MeV[ll, 12]. We find that there still remains a significant model 
dependence (Ag/t-dependence) for A ~ 400 MeV (we did not show it), and that the 
model dependence of the central-type EFT(7r)-based nuclear current almost disappears at 
A = 200 MeV (Fig. 5). Furthermore, we observe the model dependence in the tensor-type 
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Figure 5: Effective two-body operators 
(A = 200 MeV), Ol%, for the deuteron 
S-wme state from EFT(7r)-based 02,55- 
The diagonal momentum-space matrix 
elements are shown. The lower three 
dotted curves represent the original 
02,55; from the bottom, Aeft = 800, 600 
and 500 MeV. The upper three curves 
represent the 0*2 Js- The solid, dashed 
and dash-dotted curves arc the effec- 
tive two-body currents obtained from 
EFT(7r) with Kjt = 500,600 and 800 
MeV, respectively. 



Figure 6: Effective two-body operators 
(A = 200 MeV), Oafi^, for the deuteron 
D-wave state from EFT(7r)-based O^^sd- 
Three dotted curves represent the origi- 
nal 02,sD-i from the bottom, for Ae/t = 
500, 600 and 800 MeV. The other three 
curves represent the O^^lj^- The sohd, 
dashed and dash-dotted curves are for 
k^ft = 500,600 and 800 MeV, respec- 
tively. Scale of vertical axis is ten times 
as small as that of Fig. 5. The other 
features are the same as those in Fig. 5. 



operator even at A = 200 MeV (Fig. 6). This result has been expected because all of 
the operators have the same pion-range mechanism. Since the EFT(7r)-based operators 
consist of the same one-pion-exchange currant and the contact currents, they look almost 
the same only when we observe the operators with a resolution insensitive to details of 
mechanisms other than the one-pion-exchange mechanism. A resolution of A ~ 400 MeV 
is still sensitive to the short-range mechanisms, while that of A ~ 200 MeV is not. 

In Figs. 7 and 8, we plot the diagonal part of the effective two-body operators at A = 
200 MeV from the phenomenological models for the initial S- and D-wme deuteron states, 
respectively. Wc also plot the "bare" operators, 02,ss{k, k) and 02,sD{k, k) from the 
models and EFT(7r) in the figures. We find in these figures that the model dependence still 
clearly remains at this scale. This is because the starting operators are model dependent 
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Figure 7: Effective two-body operator 
^2% = 200 MeV) for the deuteron 
S-wave state from various models of 
01,88 ■ The diagonal momentum-space 
matrix elements are shown. The lower 
three dotted curves represent the bare 
02,55; from the bottom, EFT(7r), Model 
I and II. We show EFT(7r) for Kft = 600 
MeV. The dotted curve for Model II 
is almost on the zero axis. The upper 
three curves represent the Og'^^ and the 
solid, dashed and dash-dotted curves are 
the effective two-body currents obtained 
from Model I, II and EFT(7r). 



Figure 8: Effective two-body operator 
= 200 MeV) for the deuteron 
D-wave from various models of 02,5d- 
The three dotted curves represent the 
bare 02,5d; from the bottom. Model 
II, Model I and EFT(7r). The other 
three curves represent the O^^Sd- The 
solid, dashed and dash-dotted curves are 
the effective two-body currents obtained 
from Model I, Model II and EFT(7r). 
Scale of the vertical axis is ten times as 
small as that of Fig. 7. The other fea- 
tures are the same as those in Fig. 7. 



even on the one-pion exchange mechanism. 

4.2 Effective current operators for A = 70 MeV 

We now discuss our result of the WRG analysis at A = 70 MeV. In Table 1, we present 
our result for the ratio of the amplitudes from the one- and two-body operators. We 
find a significantly reduced contribution (69.22%) from the effective one-body operator at 
A = 70 MeV because the cutoff value is now rather small and thus the large portion of the 
contribution from the bare one-body operator is renormalized into the effective two-body 
operators. We also find that most of the contribution from the tensor-type two-body 
operator is renormahzed into the central-type two-body operator at A = 70 MeV^^. In 

^^In a rather small model space, the tensor nuclear force also plays an unimportant role, which is 
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Figure 9: Effective two-body axial- 
vector current operator O2ISS — 
MeV) for the initial deuteron S-wme. 
The dotted curve (almost on the zero 
line) is the same as the solid curve in 
Fig. 7. The other three curves rep- 
resent the effective currents for A = 
70 MeV. The sohd, dashed and dash- 
dotted curves are respectively derived 



Figure 10: Effective two-body axial- 
vector current operator 02^50 (-^ — 
MeV) for the initial deuteron D-wme. 
The dotted curve (almost on the zero 
line) is the same as the solid curve in 
Fig. 8. Scale of the vertical axis is ten 
times as small as that of Fig. 9. The 
other features are the same as those in 
Fig. 9. 



from Model I, Model II and EFT(7r) with 
Aeft = 600 MeV. 

Figs. 9 and 10, we show our results of the effective operators at A = 70 MeV for the 
initial deuteron S- and D-wavc states, respectively. We find again that the curves of the 
effective two-body operators suddenly change at k = A — q/2 = 55 MeV, as discussed 
earlier, due to the renormalization from the bare one-body operator. Since the three 
EFT(7r)-based operators become to be very similar to each other already at A = 200 
MeV, we show an evolution of one (Ag/t = 600 MeV) of them. With this resolution of 
the system, the dependence on modeling the details of the small scale physics is not seen 
any more. Therefore, the three different nuclear operators give essentially the same GT 
matrix elements in the kinematical region within the model space of A = 70 MeV. 

We mention here the on-shell energy dependence of the effective nuclear current oper- 
ator. As discussed in the section 2, an effective nuclear current obtained with the WRG 
equation [Eq. (2)] acquires a dependence on both the initial and final on-shell energies. 

discussed in Appendix C. 
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I ■ 1 1 

50 
t[MeV] 

Figure 11: The simulation of the effective two-body axial- vector current with A = 70 MeV 
using the contact current operators. The solid curve is the same as the solid curve in Fig. 9. 
The dotted curve simulates the solid curve using the contact current operator with no 
derivative. The dashed curve additionally has the contact current with two derivatives. 

In our case, the on-shell energy for the initial state (the deuteron) is fixed. Therefore, 
we examined the dependence of the effective two-body current on the on-shell energy of 
the final scattering state. We found a significant dependence on the on-shell energy for 
this model space (A = 70 MeV), which indicates an importance of the on-shell energy 
dependence of a NEFT-based operator defined in a small model space. 

4.3 EFT(/)-based current operator from effective current operator 

We simulate the obtained low-momentum current operator with A = 70 MeV using 
the EFT(^)-based parameterization [9, 25] 

0Ss''\k',k) = L,^ + ^ie + k-). (11) 

It is clear that the part of the effective operator after k — A — q/2 cannot be well 
simulated by this parameterization. We therefore simulate an effective operator which 
does not include the contributions from the diagrams. Figs. 4(a) and (b). The result of 
the simulation is shown in Fig. 11. Even the one-parameter fit is fairly good, and the 
two-parameter fit yields an almost perfect simulation. We note that this contact current, 
whose parameters have been fixed by fitting to the diagonal components, also simulates 
well the off-diagonal components of the effective current operator. The EFT(^)-based 
current operator can be obtained from the models or EFT(7r) in this way, and this may 



15 



be understood relation between them. 

Even EFT(^)-based current operator was determined from the models (or EFT(7r)) us- 
ing the WRG equation, it is still missing the contributions from the diagrams of Figs. 4(a) 
and (b). Because these contributions are more (less) for a larger (smaller) momentum 
transfer, validity of the operator is limited to a small momentum transfer region. For 
the ud upn reaction at E„ = 20 MeV and Tnn{= E') = 1 McV, the allowed region 
of the momentum transfer is 15 < g < 30 McV. The contributions from the diagrams of 
Figs. 4(a) and (b) to the GT matrix elements are 1.8%, 2.8% and 5.1% for q = 15, 20 and 
30 MeV, respectively. When a few percent level precision is required, this contribution is 
significant. Although there is still a kinematical region where the omission of the "jump- 
up" part is a good approximation, the EFT(|f)-based current operator constructed in the 
previous paragraph does not accurately predict the ud — > upn reaction cross sections in 
the solar neutrino energy regime {Ej, <^ 20 MeV). One useful idea to overcome the diffi- 
culty can be found in the calculations of More Effective Effective Field Theory (MEEFT) 
(or EFT*) [5], which is discussed in the next paragraph. 

MEEFT has been successfully applied to calculations of reaction rates for electroweak 
processes in few-nucleon system (see Ref. [5] and references therein), and is essentially 
equivalent to the proper NEFT as discussed in Ref. [26]. In the MEEFT calculation of 
a matrix element, the bare one-body current operator is sandwiched by wave functions 
obtained with a phenomenological nuclear force. The high-momentum components of 
the bare one-body operator have neither been integrated out nor been renormalized into 
the effective two-body operator. Meanwhile, the cutoff (Ag/t) is introduced only in 
the two-body operators. If we follow the manner of MEEFT, we would not have the 
"jump-up'' structure in the effective two-body currents. Thus we can safely parameterize 
the effective two-body operator using the NEFT-based operator and use it in a wider 
kinematical region. 

5. Conclusion 

In this work, we extended our previous renormalization group analysis of the nuclear 

force to the nuclear current o])erator. We derived the Wilsonian renormalization group 

^^A phenomenological nuclear force has a cutoff typically larger than a cutoff introduced in a NEFT- 
based nuclear force. Strictly speaking, therefore, (very) small effects of high momentum components 
of the one-body operator outside the model space should be renormalized into the effective two-body 
operator, even in the case in which we use a phenomenological nuclear force. 
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(WRG) equation for the current operator and used it for the analysis of the operators. 

In our RG analysis, we studied the five nuclear axial- vector currents, which are from the 
two phenomenological models and EFT(7r) with the three cutoff values, associated with 
the GT matrix element for the vd upn reaction in the solar neutrino energy regime. 
The original operators are different from one another due to the difference in describing 
the details of the small scale physics; even the pion range mechanism is model dependent. 
In spite of the difference, these operators give essentially the same GT matrix element 
in the low-energy region. It was found that, as the model space is reduced, the model 
dependence gradually disappears and an essentially unique effective operator is obtained 
at the small model space with A = 70 MeV. It was noted that the three EFT(7r)-based 
operators converge to a single operator in the relatively larger model space with A = 200 
MeV because they have the same pion range mechanism. As a result, the differences 
among the original operators in the reaction mechanism, such as the deuteron D-wave 
contribution in the matrix element of the two-body operator, disappear. Our RG analysis 
clearly showed the reason why the nuclear current operators of the phenomenological 
models and those of NEFT give essentially the same electroweak reaction rate. The 
reason is that as long as reactions with the kinematics included in the small model space 
are concerned, the reactions hardly probe the details of the small scale physics which 
make a difference among these nuclear operators. 

We simulated the single effective operator obtained with the WRG equation using 
the EFT(|f)-based parameterization. A good simulation indicates that if we observe the 
operators of EFT(7r) or phenomenological models roughly enough, they look like that of 
EFT(/), not only in the strong sector but also in the electroweak sector^^; this is the 
relation among these operators. Furthermore, we can regard the EFT(]f)-bascd current 
operator as model-independent because it simulates well the unique, model-independent 
effective current operator defined in a certain small model space. However, the accurate 
simulation can be done only if we omit the "jump-up" part of the effective operators. This 
omission limits ourselves to a rather small kinematical region (small momentum transfer) 
where the omission gives negligible effects to the GT matrix element. In this regard, we 
discussed the advantage of using MEEFT to avoid the problem. 

Finally, we make some remarks here. Although we performed the RG analysis of the 

^■^If the pion range mechanism of the phenomenological nuclear current models were the same as that 
of EFT(7r), we could have also connected them through WRG, as we did for the nuclear force in Ref. [11]. 
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nuclear current operator for the ud upn reaction, the result we obtained should be 
the case for other electroweak processes in two-nucleon system such as the np d'j 
reaction [27, 28] and the pp de^v^ reaction [20, 22, 25, 29, 30]. That is, in nuclear 
current operators for those two-nucleon processes, one can also find the relation between 
models and NEFT using the WRG method. The operators for the vd vpn reaction are 
studied as an example to study the general relation. With the general result, we mention 
an impact of this work on non-nuclear physics communities, such as neutrino physics and 
astrophysics. Theoretical cross sections of various reactions in two-nucleon system have 
been often used by people in these communities as input parameters in their analysis of 
their own problems. Although the agreement in the theoretical cross sections between 
the different approaches (models, NEFT) has been a good basis for those people to rely 
on the theoretical values, the situation is not totally satisfactory because the reason for 
the agreement has not been clearly shown. Now we can provide the communities with a 
firmer basis, i.e., the equivalence relation between phenomenological models and NEFT 
from the viewpoint of RG, which we have demonstrated in this work. 
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Appendix A: Derivation of Wilsonian renormalization equation for nuclecir 
current operator 

We start with a matrix element, evaluated in the momentum space, of an operator O 
defined in a model space: 

iO) =1^1-^ Oi'^^\k', k;p',p) 4-\k) , (12) 

where ijj^°'\k) {'<p^^\k')) is the radial part of the wave function for the initial (final) two- 
nucleon state, and a {pi) specifies the partial wave. The quantities k and k' {p and p') 
are the magnitudes of the (on-shell) relative momenta of the two-nucleon system. The 
cutoff A, which is the maximum magnitude of the relative momenta, specifies the size 
of the model space for the two-nucleon states. The radial part of the current operator 
between the 13 and a partial waves is denoted by 0^f'°'\k', k;p',p), which may depend on 
the on-shell momenta and the cutoff. 
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We differentiate the both sides of Eq. (12) with respect to A, and obtain 

^ /;« ,^..(.0 ^^""(^■^^^-■^) . (13) 

Here we write the wave function as 

= ..^^-^ , ^-L— £t^vtH,. <>(.') , (14, 

where = /M with M being the nucleon mass. The quantity V^\k, k';p) is the A^A^- 
potential defined in the model space specified by A; the A-dependence of vl°'\k, k';p) is 
controlled by the WRG equation for A"A^ potential [11, 13]. It is noted, in the case of 
the deuteron wave function, that the first term in the r.h.s. of Eq. (14) does not exist, 
and that the normahzation is such that the amplitudes are the same as those of the wave 
function from Va=oo whose squared integral is normalized to be unity. Using Eq. (14), we 
can rewrite Eq. (13) as 

- /;*r^^*'r(*ovr(*',A;.')^or(A,..',.)<'(^) 
+ /;« oofjM 

Now we impose a condition that the matrix element is invariant with respect to changes 
in A, i.e., 

d{0) 



dA 



= , (16) 



then we find from Eq. (15) the Wilsonian renormalization group equation for the operator 
O: 

dOf''\k',k-p',p) M A^ , , . 

dK " 2^ ^ ' '^ \^-p'^ ^ \Kk-p,p) 

M 

+ i^,Of'^\k'.^:P\p)j^,vlr\AMp). (17) 
19 



Note that Eq. (17) is a sufficient condition of Eq. (16), but not a necessary condition. 
Nevertheless, Eq. (17) is physically appealing because this equation can also be derived 
with the path integral method, as stated in footnote 7, without concerning about whether 
Eq. (17) is a necessary condition or not. 

Appendix B: Conservation law for effective current operator 

It is interesting to study a conservation law for an effective current operator obtained 
with the WRG equation [Eq. (2)]. Suppose that we start with a model in which the 
(vector or axial-vector) current operator J'^ satisfies an equation 



(18) 



where S denotes a source term and is the momentum transfer from the external current 
to a nuclear system. (We will use as the momentum transfer in the following.) We also 
suppose that this model is defined in a model space whose size is specified by a cutoff A. 
Now we reduce the model space of the current operator using the WRG equation up to a 
cutoff A. We obtain 

where J^^^ is the effective current operator defined in the model space with A. The full 
Hamiltonian and the iViV-interaction are respectively denoted by H and V ^ both of which 
are defined in the model space with A. The quantities E and E' are the kinetic energies 
of the two-nucleon system in the initial and final states, respectively. The projection 
operators 77 and A are defined by 



X 



J (2^ 



k)(k 



k 



< A 



r (Pk 
J (2^ 



k)(k 



A < 



k 



<A, 



(20) 
(21) 



where \k) represents the free two-nucleon states with the relative momentum k. 

Now, with the Eqs. (18) and (19), the equation which the effective current satisfies is 



ff 



^eff 



with 



Seff -v{s + SX^XV + VX^XS + ^A^A^A^A^) , . 



(22) 



(23) 
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Therefore, if the current in the starting model is conserved (S = in Eq. (18)), then 
the effective current J^^j is also conserved {Seff = in Eq. (22)). As has been done in 
the text, we simulated the effective current with the use of the contact currents. In such a 
case, the contact currents violate Eq. (22) to a certain extent. However, as we have seen, 
the simulation is quite accurate, up to the "jump-up" part, and therefore the violation of 
the equation is very small. 

Appendix C: Tensor nucleeir force in pionless effective field theory 

In describing low-energy two-nucleon system, two types of nuclear effective field the- 
ories (NEFTs) have been often used. One of them considers the nucleon and the pion 
exphcitly [EFT(7r)] while the other treats only the nucleon as an explicit degree of freedom 
[EFT(/)]. Although both NEFTs successfully describe low-energy NN system, there is 
a distinct difference between them in dealing with the tensor force. In EFT(7r) with a 
cutoff regularization scheme, the tensor force which comes from the one-pion-exchange 
potential (OPEP) is considered to be a leading order (LO) term, and is iterated to all 
orders by solving the Schrodinger equation. Thus, the tensor force in EFT(7r) is LO and 
treated non-perturbatively^^. In contrast, treatments of the tensor force in EFT(/) are as 
follows. In the cutoff scheme, the tensor force is a ncxt-to-leading order (NLO) term and is 
treated non-perturbatively (iteration to all orders). In the KSW counting scheme[31], the 
tensor force is a next-to-next-to-leading order (N^LO) term and is perturbatively treated. 
Thus, the tensor force is considered to be less important in EFT(/), which is based on 
the counting rules used. However, it is still interesting to examine whether the tensor 
force indeed has different importance in EFT(7r) and EFT{f). This examination can be 
done quantitatively by deriving VEFTifi from VEFT{n), for which the renormalization group 
plays an essential role {VEFT{n) (^FT(;t)) is the nuclear force based on EFT(7r) {EFT{f))). 
This is what we discuss in this appendix, using a result given in our previous paper [11]. 

We use Veft{tt) and VEFT{f) presented in Ref. [11] in the analysis here. More specifically, 

we use V^j^r(^) which consists of the OPEP plus contact interactions with zero, two and 

four derivatives , and V^|'T(;f) which consists of contact interactions with zero, two and four 

derivatives. The cutoff values, which specify the size of the model space for the nuclear 

forces, are A = 200 MeV for ^i^^M ^nd A = 70 MeV for V^^^j^^^y In Ref. [11], the 

parameterization of these nuclear forces is given in Eq. (2.14) and Appendix B, and the 

Although a perturbative treatment of the OPEP was proposed in the KSW counting scheme[31], it 
turned out that the OPEP has to be treated non-perturbatively [32]. 
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numerical values of the parameters are given in Table I. The two interactions, V^|^7-(^) and 
'^EFT{f) > related to each other through the Wilsonian renormalization group equation 
(Eq. (2.3) of the reference). 

Now we focus our attention on the ^Si-^Di partial wave, because we are interested in 
the importance of the tensor force. We calculate the NN scattering phase shifts and the 
deuteron binding energy using V^^pj,(^^^^ and V^Fr(^)- In order to study the importance of 
the tensor force, we switch on and off the tensor force included in Vj^pj.^^-^ and V^^t(^)- 
The result is shown in Table 2. We can observe that the tensor force in Vpprp^^^^ gives 
only a tiny effect to the observables. On the other hand, the tensor force in Vppj,^^^ plays 
an unnegligible role, as expected. In this way, we conclude that the tensor force has 
different importance in EFT(7r) and EFT(^), which supports the counting rules used so 
far. 
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A 


Model 


(IB) 


(2B) 


(2B) 


(MeV) 




S + D (%) 


S + D (%) 


S{%) 


Di%) 




Model I 


98.87 


1.13 


-0.33 


1.46 




Model II 


98.87 


0.72 


0.00 


0.72 


oo 


Aeft = 500 MeV 


98.87 


1.42 


-1.18 


2.60 




Ae/t = 600 MeV 


98.87 


1.41 


-2.00 


3.41 




Aeft = 800 MeV 


98.87 


1.40 


-3.26 


4.65 




Model I 


97.05 


2.96 


2.85 


0.12 




Model II 


97.05 


2.56 


2.50 


0.06 


200 


Aeft = 500 MeV 


97.05 


3.25 


3.05 


0.20 




Aeft = 600 MeV 


97.05 


3.25 


3.04 


0.21 




Aeft = 800 MeV 


97.05 


3.23 


3.01 


0.22 




Model I 


69.22 


30.79 


30.77 


0.03 




Model II 


69.22 


30.39 


30.36 


0.03 


70 


Aeft = 500 MeV 


69.22 


31.08 


31.05 


0.03 




Aeft = 600 MeV 


69.22 


31.08 


31.05 


0.03 




Aeft = 800 MeV 


69.22 


31.06 


31.04 


0.03 


: Contributions (ratio) 


to the Gamow-Teller 


(GT) 


matrix 



{^So\0\^Si{^Di)), from each component of the current operators. The model space reduc- 
tion using the WRG equation introduces the sharp cutoff and the values of it are shown 
in the first column. Each current operator is obtained from Model I and II, and EFT(7r) 
with Aeft = 500, 600, 800 MeV, shown in the second column. The third and fourth 
columns give the contribution from the one-body (IB) and the two-body (2B) operators, 
respectively. The 2B contributions are from the deutcron '^5*1 and ^Di waves, which are 
separately given in the fifth and sixth columns (denoted by S and D). For more details, 
see the text. 
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Table 2: The importance of the tensor force in NN observables. The first column indicates 
the nuclear force used, while the second one specifies whether the tensor force contained in 
the nuclear force is turned on or off. The third (fourth) column presents the phase shifts 
of the A^A^ '''Si-'^'Di partial wave scattering for p = 10(30) MeV; p is the on-shell relative 
momentum of the two nucleons. The fifth column gives the deuteron binding energy. 







tensor force 


"p=10MeV 


"p=30MeV 


B.E. 




^ EFT(-k) 


ON 


164.47 


137.19 


2.225 




= 200 MeV) 


OFF 


163.46 


135.15 


1.835 




^EFT(f) 


ON 


164.47 


137.35 


2.223 


(A 


= 70 MeV) 


OFF 


164.46 


137.34 


2.217 
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